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Abstract — In this work, a non linear two-tank interacting 
level process is taken-up for study. The mathematical model of 
the two-tank interacting process is derived. A conventional PI 
controller is designed based on process reaction curve with ZN 
technique. To improve the performance of PI controller, a 
Particle Swarm Optimization (PSO) based PI controller PSO-PI 
is proposed and implemented. The servo and regulatory 
responses with PI and PSO-PI are obtained and discussed. From 
the results and performance measures, it is observed that 
PSO-PI is producing better results than conventional PI 
controller. 

Index Terms — Two- tank interacting process, PI controller 
and PSO-PI. 

I. INTRODUCTION 

The process industries require the liquids to be pumped, 
stored in tanks and then pumped to another tank. Many times 
the liquid will be processed by chemical or mixing treatment 
in the tanks, but always the level of the fluid in the tanks must 
be controlled. A level that is too high may upset the reaction 
equilibrium, cause damage to the equipment or result spillage 
of valuable or hazardous material [1]. If the level is too low, it 
may have bad consequences for the sequential operations. 
Hence, control of liquid level is an important and common in 
process industries. Conventional controllers are widely used 
in industries since they are simple robust and familiar to the 
field of operator. The most basic and pervasive controller 
algorithm used in the feedback control is proportional integral 
controller algorithm. PI controller is widely used control 
strategy to control most of industrial automation process 
because of its remarkable efficiency and simplicity. Hence, a 
mathematical model is derived and a simulation is carried out 
for the given mathematical equation. The conventional PI 
controller parameters are designed based on 
Ziegler-Nicholas. However, it will not give satisfactory 
response for change in operating point. Hence, the Particle 
Swarm Optimization technique is proposed and implemented 
to optimize PI controller parameters [4]. In this work, a 
simple performance criterion in time domain is proposed for 
evaluating the performance of a PSO-PI controller that is 
applied to a interacting two-tank process. 

II. PROCESS DESCRIPTION 


valves, six process tanks, two overhead tanks, two differential 
pressure transmitters, five level transmitters and rotameters 
.Instrumentation panel consists of two PID controllers, main 
power supply switch, pump switches, motorized control valve 
switches and auxiliary switches for individual components. 

III. Working principles 

Fluid level in the tank is measured by level transmitter (LT). 
Output of LT is given to the data acquisition setup. It consists 
of ADC and DAC. The differential pressure level transmitter 
(DPLT) measures the flow by sensing the difference in level 
between the tanks. The DPLT then transmits a current signal 
(4-20mA) to the I/V converter. The output of the I/V 
converter is given to the interfacing hardware associated with 
the personal computer (PC). Control algorithms are 
implemented in Lab view software. It compares and takes 
corrective action on the motorized control valve. Based on the 
valve opening flow rate is manipulated. Rotameter can 
visualize the flow rate. 



Fig.l. Experimental setup of a two-tank interacting process. 


The controller compares the controlled variable against set 
point and generates manipulated variable as current signal 
(4-20mA) [2]. Here the controlled variable is the level (h 2 ) 
and the manipulated variable is the flow rate 
(q in ). Specification dimension are tabulated in Table I and 
Table II respectively. 


Fig.l shows the photograph of the laboratory level process 
station. It consists of three pumps, two motorized control 
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TABLE I 

SPECIFICATION OF TWO-TANK INTERACTING 

PROCESS 


COMPONENTS 


SPECIFICATIONS 
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MOTORISED 

CONTROL 

VALVE 

Flow rate 
Characteristics 
Valve action 

50LPH 

Equal% 

motorized 

control 

ROTAMETER 

Type 

Range 

float material 

variable area 
(O-IOO)LPH 
i.e (0-1666 
3 /min) 

SS 316 

PUMP 

RPM 

Discharge 

Voltage 

4500 

1000(LPH) 
220/230 volts 
AC&DC 

PROCESS TANK 

Capacity 

Height 

Diameter 

3 litres 
300 mm 
120 mm 

LEVEL 

TRANSMITTER 

Input 

Height 

Type 

24V DC 
0-400 mm WC 
F capacitance 

DIFFERENTIAL 

PRESSURE 

TRANSMITTER 

Input 

Supply 

Output 

0-400mm H 2 0 
24V DC 
4-20mA at 24V 
DC 


TABLE II 

DIMENSIONS AND VARIABLES FOR TANK1 AND 

TANK 2. 


Parameter 

Dimension 

Area(Ai,A 2 ) 

113.0973 cm 2 

Height(max) 

25 cm 

Diameter 

12 cm 

Inflow rate(MV),q in 

0-1666 cm /min 

Process Variable (PV),h 2 

0-25 cm 

Flow rate (LV),q L1> q L2 

0-500 cm 3 /min 

G 

Gravity (9.81N/m 2 ) 

aiand a 2 

Area of the pipe outlet 
(ai=3. 5735cm 2 , 
a 2 =3. 9012cm 2 ) 

c d 

Discharge co-efficient 
( c d =0.08) 


tank2 (cm), respectively Both tanks have the same cross 

2 2 

sectional area A^cm ) and A 2 (cm ), q L1 is the inflow of tankl 
as load disturbance (cm /min) and q L2 is the inflow of tank2 
as load disturbance (cm /min) [land 2]. 

qin(mv) qu (Load) qL 2 (Load) 



tankl tank2 

Fig.2 .Two-tank interacting process 


Mass balance equation for tank 1 

dh, 

A i at = q > » ~ ? > 

Assuming non -linear resistance to flow : 
1 

Mass balance equation for tank 2 



-q 0 


Assuming non-linear resistance to flow: 



y 2 

A_ 







C _cf 1 

A 


MV- 


C 

d 2 

A 


2g/t, 



MV- manipulated variable, PV-process variable LV- 
,6. Load variable 
inti 

Fig.2 shows the process consists of two interacting liquid 
tanks. The volumetric flow into the tankl is q in (cm Ymin) , the 
volumetric flow rate from tankl to tank2 is qi(cm /min), and 
the volumetric flow rate from tank2 is q 0 (cm /min) [2] . The 
height of the liquid level is hi (cm) in tankl and h 2 (cm) in 


IV. CONTROLLER TUNING 

A. Process reaction curve for two -tank interacting process 

Fig. 3 shows the simulated open loop response of 
interacting process. The level (h 2 ) changes from 0 to 8.5 cm, 
when applying a step input in q in (50* 16.66cm Ymin) also the 
level (hd changes from 0 to 18 cm due to interaction. The 
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simulated process reaction curve(PRC) of h2 for step change 
in q in for ±499. 8cm /min is shown in Fig.4. 



Fig.3. Simulated open loop response of hi and h 2 of 
interacting Process. 



Time(sec) 


Fig.4. Simulated PRC of h 2 for step change in q in for 
+499.8cm 3 /min. 

From the fig.4 the transfer functions are obtained and 
tabulated in Table III. From the average transfer function, the 
controller parameters are obtained using Z-N tuning rule [2] . 
For two-tank interacting process the PI controller parameters 
are tabulated in Table IV. 


TABLE III 

TRANSFER FUNCTION MODEL OF TWO-TANK 
INTERACTING PROCESS 


Step 

Input q in 

Transfer Function 

Average Transfer 
Function 

Positive Step 
Input, q in 

0.81578 

„- 17 j 

0.65049 

935 + 1 “ 

Negative Step 
Input, q in 

0.4352 

p-l-TS 

97.5s + 1 

1025 + 1 " 
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INTERACTING PROCESS USING Z-N METHOD. 


Mode 

K c 

r t (jec) 

PI 

7.9351 

56.61 


C. Particle swarm optimization (pso) algorithm 

Particle swarm optimization (PSO) algorithm, which is 
used for optimizing the PSO difficult numerical functions and 
based on metaphor of human social interaction, is capable of 
mimicking the ability of human societies to process 
knowledge [4]. It has roots in two main component 
methodologies: artificial life (such as bird flocking, fish 
schooling and swarming) and evolutionary computation. 
Although the PSO algorithm is initially developed as a tool 
for modelling social behaviour, it has been applied in 
different areas. 



Fig.5 . PSO based PI controller to control the level of 


two- tank interacting process. 


The block diagram of PSO based PI controller is shown in 
Fig 5. Also the flow chart to optimize the PI controller is 
shown in Fig 6. 



Fig 6.Flow chart PSO-PI 


A .Z-N tuning technique 

Based on the average transfer functionthe value of k c and 
f are calculated using Z-Ntechniques.The values of kc and f 
are tabulated in Table IV. 

TABLE IV 

PI CONTROLLER SETTINGS FOR TWO-TANK 


The equation for position and velocity 

V n+ i = v n +Cirandl( )*(p bes t, n -current position)+ c 2 randl( 
)*(gbest,n-current position) 

X (t +1) = x ld (t) + v ld (t +1) 


(5) 

( 6 ) 
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D. The PSO learning algorithm is represented by 
following steps 

Stepl Select the number of iterations (n) and the PSO 
learning rate (ci, c 2 ). 

Step2 Randomly generate position vector X and 
associated velocity V of all particles in the 
population. 

Step3 Calculate each individual’s fitness value, and 
then compare each individual’s evaluation value 
with best global particle value g bes t the personal 
best value p best . Finally select the new best 
value (g b 

est and p b e S t) 

Step4 For every particle, update its own velocity and 
position value. 

Step5 If iterations, then go to exit, otherwise go to 
step 3. 

Step6 The best particle’s value will be selected as the 
finial parameter set to form desired PI 
controller. 


During simulation the values for PSO algorithm 
considered are tabulated in Table V. 


TABLE V 

PARAMETER OF PSO ALGORITHM 


Parameter 

Values 

Number of Iteration 

20 

Dimension 

3 

Swarm Size 

90 

Correction Factor 

Cl 

1.2 

C 2 

1.4 


The optimized PS O-PI parameter using equation 5 and 6 
are tabulated 

TABLE VI 


PARAMETERS OF PI AND PSO-PI CONTROLLERS 


Mode 

K c 

T^sec) 

PI 

7.9351 

56.61 

PSO-PI 

2.03791 

0.1005 


V. SIMULATION RESULTS 

A. Servo responses of level h 2 with PI and PSO-PI 
controllers. 


The servo responses for two tank level interacting process 
for step change in h 2 from 6 to 8 cm, from 8 to 6cm is shown 
in fig. 7 by implementing PI and PSO-PI. The corresponding 
change in h x and inlet flowrate q in are also shown in fig. 8 and 
fig. 9, respectively. The performance measures for step change 
in h 2 from 6 to 8 cm is tabulated in Table VII. From the results 
and Table VII, it is observed that the PSO-PI gives better 


performances than PI controller. 



Fig.7. Servo response of h 2 with PI and PSO-PI 
controllers. 




Fig 9. Response of PI and PSO-PI output qin for 
servo response. 
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TABLE VII 

COMPARISON OF PERFORMANCE MEASURES 
LEVEL WITH PI AND PSO-PI FOR SERVO 
RESPONSE (6-8) CM 


Servo response for h 2 (6-8)cm 


t s (sec) 

% overshoot 

ISE 

PI 

1988 

12.825 

353.4 

PSO-PI 

978 

33.66 

254.2 


B. Regulatory response for Positive step change load 
variable qLi(+) 

A sudden load disturbance of +10% is given in inlet 
flow rate of tankl at 4000 sample from q L1 (+) . Due to this 
level in tq increases for PI from 13.3 to 14.5 cm and PSO it 
increases from 13.3 to 15.4 as shown in fig.ll. The level h 2 
also increases for PI from 6 to 6.5 cm and PSO it increases 
from 6 to 6.9cm as shown in fig 10. The PI and PSO-PI 
Controller takes necessary action to reduce the flowrate q in 
from 420 to 350 (ref fig 12) In order bring back the level h 2 
back to 6cm as shown in fig 10. 




Fig 10. Regulatory response of h 2 with PI and PSO due 
to load variation in +10% from q L1 . 


CO 


O 

c 

C7^ 



4000 


4500 5000 

time (sec) 


5500 


Fig 12. Response of PI and PSO output(q in ) for load 
variation in +10% from q L1 


C. Regulatory response for negative step change load 
variable qLi(-) 


A sudden load disturbance of -10% is given in inlet flow 
rate of tankl at 4000 sample from qu(-) . Due to this level in 
hy decreases for PI from 13.3 to 12.1 cm and PSO it 
decreases from 13.3 to 11.1 as shown in fig 14. The level h 2 
also decrease for PI from 6 to 5.5 cm and PSO it decreases 
from 6 to 5.1 cm as shown in fig 13. The PI and PSO-PI 
Controller takes necessary action to increase the flowrate q in 
from 425 to 500 (ref fig 15) in order bring back the level h 2 
back to 6cm as shown in fig 13. 




4000 4500 5000 5500 


time(sec) 



Fig 11. Regulatory response of fq with PI and PSO-PI due to 
Load variation in +10% from q L1 


Fig 13. Regulatory response of h 2 with PI and PSO-PI due to 
load variation in -10% fromq L1 



4000 4500 5000 5500 


time(sec) 

Fig 14. Regulatory response of hi with PI and PSO-PI due 
to load variation in -10% from q L1 
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Fig 15. Response of PI and PSO-PI output (q in ) for load 
variation in -10% from q L1 

TABLE VIII 

COMPARISON OF PERFORMANCE MEASURE OF 
LEVEL WITH PI AND PSO FOR REGULATORY 

RESPONSE. 


Regulatory Response h 2 (6-8)cm 

Controller 

+10% from q L1 

-10% from q L1 

t s (sec) 

ISE 

t s (sec) 

ISE 

PSO-PI 

1032 

2822 

1046 

2821 

PI 

1460 

9316 

1275 

9316 


D. Regulatory response for positive step change load 
variable qL 2 (+) 


A sudden load disturbance of +10% is given in inlet flow 
rate of tankl at 4000 sample from qL 2 (+) • Due to this level in 
tq increases for PI from 13.3 to 13.5 cm and PSO it increases 
from 13.3 to 13.6cm as shown in fig 17. The level h 2 also 
increases for PI from 6 to 7 cm and PSO it increases from 6 to 
6.9cm as shown in fig 16. The PI and PSO-PI Controller 
takes necessary action to reduce the flowrate q in from 420 to 
200 (ref figl 8) in order bring back the level h 2 back to 6cm as 
shown in fig 16. 




time(sec) 

Fig 16. Regulatory response of h 2 with PI and PSO-PI 
due to load variation in +10% from 


4L2 



. Fig 17. Regulatory response of fq with PI and PSO-PI 


due to load variation in +10% from q L2 



time(sec) 

Fig 18. Regulatory response of tq with PI and PSO-PI 
due to load variation in +10% from q L2 


E. Regulatory response for Negative step change 
load variable qL 2 (-) 

A sudden load disturbance of -10% is given in inlet flow 
rate of tankl at 4000 sample from qL 2 (-) • Due to this level in 
tq increases for PI from 13.3 to 15.6 cm and PSO it increases 
from 13.3 to 15 as shown in fig 20. The level h 2 also decreases 
for PI from 6 to 5.6 cm and PSO it decreases from 6 to 5.8cm 
as shown in fig 19. The PI and PSO-PI Controller takes 
necessary action to increase the flowrate q in from 420 to 750 
(ref fig21) in order bring back the level h 2 back to 6cm as 
shown in fig 19. 




Fig 19. Regulatory response of h 2 with PI and PSO-PI due 
to load variation in -10% from q L2 
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Fig 20. Regulatory response of hi with PI and PSO-PI due 
to load variation in -10% from q L2 



Fig 21. Regulatory response of q in with PI and PSO-PI 
due to load variation in -10% from q L2 . 

TABLE IX 

COMPARISON OF PERFORMANCE MEASURE OF 
WITH PI AND PSO-PI FOR REGULATORY 
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Regulatory Response h 2 (6-8)cm 

Controller 

+10% from q L2 

-10% from q L2 

ts(sec) 

ISE 

ts(sec) 

ISE 

PSO-PI 

1036 

2824 

1045 

2826 

PI 

1460 

9321 

1277 

9311 


VI. CONCLUSION 

The non-linear mathematical model of a two-tank 
interacting process was derived. To improve the performance 
of closed loop control a PSO-PI was designed and 
implemented for a two-tank interacting process. The servo 
and regulatory responses were obtained with PSO-PI. The 
performances of PSO-PI were compared with that of 
conventional PI controller in simulation. The performance 
measures were tabulated. It is observed that the PSO-PI gives 
is better performance on terms of less integral square error, 
faster settling time and without oscillation. 
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